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Abstract. We present a model for estimating the blazar contribution 
to the extragalactic 7-ray background (EGRB) in the EGRET energy 
range. It is based on inverse-Compton scattering as the dominant 7-ray 
production process in the jets of flat spectrum radio quasars (FSRQs) 
and BL Lac objects, and on the unification scheme of radio-loud AGN. 
Depending on the evolutionary behaviour of radio-loud AGN, our esti- 
mate for the point source contribution to the extragalactic background 
emission ranges from about ~ 20% — 40% for a redshift drop off at 
Zmax = 3 to around ~ 35% — 85% for a cutoff at Z max = 5. Roughly 70%- 
90% of the AGN contribution to the EGRB would be due to BL Lacs/FR I 
objects. However, in case of BL Lacs our model parameters are not en- 
tirely in agreement with the results from studies at radio wavelengths 
indicating systematical uncertainties whereas our estimate for the FSRQ 
contribution seems to be consistent. Thus we end up with two possibilities 
depending on whether we underpredict or overpredict the BL Lac contri- 
bution: either unresolved AGN can not account for the entire EGRB or 
unresolved BL Lacs produce the observed background. 



Introduction 



One of the most important results of the Energetic Gamma Ray Experiment 
Telescope (EGRET) in the field of extragalactic astronomy is the measurement 
of an isotropic diffuse 7-ray background over three decades of photon energy 
(Sreekumar et al. 1998). The differential spectrum can be described by a power 
law with photon spectral index a = 2.10 ± 0.03. The intensity of the EGRB 
integrated above 100 MeV is I = (1.45 ± 0.05) • 10~ 5 photons cm 2 s 1 sr . 
In this paper we follow a new ansatz for modeling the blazar contribution to 
the EGRB, which is based on the unification scheme of radio- loud AGN and a 
particular 7-ray emission model, namely external inverse-Compton scattering of 
accretion disk photons (Dermer & Schlickeiser 1993). 
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2. The model 



According to the unification scheme of radio-loud AGN blazars are the beamed 
sub-class of Fanaroff-Rilcy (FR) radio galaxies (i.e. Urry & Padovani 1995) 
with the most probable parent population of FSRQs being the FanarofF-Riley 
II (=FR II) radio galaxies and BL Lacs being the beamed population of FR I 
radio galaxies (i.e. Padovani & Urry 1992; Urry, Padovani & Stickel 1991). For 
a homogeneous distribution of viewing angles and an isotropic distribution of 
these objects throughout the universe, the 7-ray emission from radio-loud AGN 
observed under a large viewing angle and/or at a large distance is expected to 
fall below the EGRET-sensitivity limit. The integration of the 7-ray intensity 
emitted by those objects yields the point-source contribution to the EGRB. 
We base the calculation of the non-thermal 7-ray emission of blazars on the 
leptonic model of Dermer & Schlickeiser (1993) and extend it by considering 
not only instantaneous but also continuous electron injection. In particular 
we consider a power law with index s=2 for the injected electron spectrum 
and a fixed luminosity for the Shakura-Sunyaev-accretion disk luminosity for all 
sources. The calculation of the template source spectrum is described in detail 
in Miicke & Pohl (1998). 

Using a power law distribution for the energy injected into the jets of AGN, 
PFR(Ei n ) = 10 4?) crg ( 10 4iie rg )~ <5 (E in > E in>1 ), we calculate the time and energy 
integrated flux (from 100 MeV to 10 GeV) from BL Lac/FR I and FSRQ/FR II 
objects. The low energy cut-off Ej n 1 and the normalization po of the distribution 
PFn(Ein) are the free parameters of the model. 

We use a Monte-Carlo method to model the Log N-Log S and the redshift 
distribution of resolved BL Lacs and FSRQs which are our model constraints. 
Objects which are brighter than the sensitivity limit determine the resolved 
Log N-Log S, while sources with flux values below that limit contribute to the 
EGRB. The unresolved blazar contribution to the EGRB then follows by sum- 
ming up all unresolved AGN per solid angle of the whole survey. Our calculations 
consider a redshift drop off at Z max = 3 and Z max = 5. Note that the 7-ray 
emission in the EGRET energy range is moderated by pair production on the 
IR-background if it originates at redshifts Z > 3 (Stecker et al. 1992). The 
observed power law spectrum of the EGRB up to 100 GeV therefore excludes 
significant contributions from sources at redshifts of Z > 3. 

3. Results and conclusions 

The Log N-Log S and redshift distributions of both, FSRQs and BL Lacs, were 
compared with our model distributions by using a Kolmogorov-Smirnov test. A 
compilation of the derived values for the free parameters can be found in Table 
1 for a redshift cutoff at Z max = 3 and in Miicke & Pohl (1998) for Z max = 5. 
With the knowledge of Ej n 1 and po the distributions of the injected energy of 
BL Lacs/FR I and FSRQs/FR II can be constructed. In general FR I/BL Lacs 
are more numerous than FSRQ/FR II, but less luminous. The total luminosity of 
all FR I/BL Lacs is similar to that of all FSRQ/FR II though. As a consequence 
the background contribution of unresolved AGN is dominated by FR I/BL Lacs 
(~ 70% - 80%), while resolved radio-loud AGN are mainly FSRQs. With our 
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model we can only explain ~ 20% — 40% of the observed extragalactic 7-ray 
background as caused by unresolved radio-loud AGN. If the redshift cut-off is 
assumed to occur at Z max = 5 this contribution increases to ~ 35% — 85% of 
which ~ 80% — 90% would be BL Lacs/FR I. The typical unresolved source 
in our model is either weakly beamed or has a weak intrinsic luminosity of the 
blazar jet. 

In Figure 1 we show the predicted differential power distribution S 2 dN/dS of 
FSRQs/FR II and BL Lacs/FR I for the case of instantaneous electron injection 
and a source density cut-off at Z max = 3. It indicates that the main contribution 
to the EGRB stems from sources with a flux around S ~ 10~ 8 ph cm~ 2 s" 1 . 
Our model predicts a significant flattening of the Log N-Log S distribution in 
the next two decades of source flux below the EGRET threshold and number 
counts of roughly ~ 0.1 deg -2 . 

Comparing the total number M of BL Lacs/FR I and FSRQs/FR II from our 
model fit with the observations in the radio band (see Padovani & Urry (1992), 
Urry, Padovani & Stickel (1991)) reveals a significant underestimate of the to- 
tal number of BL Lac/FR I-objects while for the FSRQ/FR II population our 
model is in agreement with the radio observations. We note that the 7-ray 
emission model used in our study may not be ideally suited for BL Lacs anyway, 
since it is based on the assumption that the accretion disk provides the target 
photons for the dominant inverse-Compton process (EIC). However, BL Lacs 
are thought to have rather weak accretion disks which may imply that the SSC 
mechanism is the dominant 7-ray radiation process. We therefore conclude that, 
depending on whether we underpredict or overpredict the BL Lac contribution, 
either unresolved A GN can not account for the entire EGRB or unresolved BL 
Lacs/FR I produce the observed background. The detection of roughly 10 2 — 10 3 
more blazars in the energy range of ~ lOOMeV - lOGeV with a telescope ~ 10 
times more sensitive compared to EGRET is expected. 
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FSRQs 


BL Lacs 




instan.inj. 


cont.inj. 


inst.inj. 


cont.inj. 


Ei„,i 


10 44.4 erg 


10 44 ' 4 erg 


10 4J.» crg 


10 4J.» crg 


PO, 


(5.9 ± 1.0) ■ 10 6 


(4.9 ±0.8) ■ 10 6 


80 ± 27 


46 ± 16 


N 


(9.43 ± 1.53) • 10 4 


(7.82 ± 1.26) • 10 4 


(2.71 ±0.90) • 10 5 


(1.56 ±0.52) ■ 10 5 


N[%\ 


26± 5 3 


33+ 3 


74+ 3 


67l 3 5 


Pns 


0.61 


0.37 


0.40 


0.36 


Pz 


0.65 


0.71 


0.96 


0.78 


Sunres 


0.90 ± 0.15 


0.99 ±0.16 


3.24 ± 1.08 


2.78 ± 0.93 


Sunres [%] 


6± 1 


7± 1 


22 ±7 


19+6 


Sunres 
Sunres , tot 


22± 4 


26+J 


78^ 


74 +3 



Table 1. The parameters of our model for Z max — 3 together with 
the KS-probability for the statistical identity between the data and 
model Log N-Log S distributions (Pns) and redshift distributions (Pz)- 
Ej ni i and po are the low energy cut-off and the normalization of the 
injected energy distribution, respectively, and are the fit parameters of 
our model. M is the total number of objects of each object class in the 
whole universe and po (in Mpc~ 3 erg" 1 ) the normalization of the distri- 
bution p(Ei n ). S unres is the intensity (given in 10~ 6 photons cm" 2 s" 1 
sr" 1 above 100 MeV and as fraction of the observed background radi- 
ation) due to the sum of unresolved BL Lacs/FR I and FSRQs/FR II 
while Sunres gives the fractional contribution of each object class. 

Sunres, tot 

The given uncertainties are derived by varying the number of resolved 
objects by ±y/N. We predict that ~ 20% - 40% of the EGRB is due 
to unresolved AGN of which ~ 70% - 80% would be BL Lacs/FR I. 
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Fig.l: Predicted Log (S 2 
(dN/dS))-Log S-distribution 
of our background model 
for BL Lacs/FR I and FS- 
RQs/FR II for the case of 
a finite Ej n -distribution and 
instantaneous electron injec- 
tion assuming a redshift cut 
off at Z = 3. The statisti- 
cal uncertainty in each bin 
lies at about 10%. The main 
contribution to the EGRB 
comes from sources at S ~ 
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